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An  air  flow  facility  was  used  to  simulate  the  flow  through  a set  of  wicket 
gates  and  the  draft  tube  typical  of  a hydroelectric  pump-turbine.  A method 
of  preventing  draft  tube  surge  was  studied,  which  consists  of  injecting  some 
high  energy  fluid  from  the  spiral  casing  into  the  draft  tube  counter  to  the 
direction  of  the  swirl.  The  results  of  these  tests  provide  design  criteria 
for  selecting  injection  nozzle  configuration  and  estimating  the  required 
injection  flow  for  a turbine  of  given  performance  characteristics.  A potential 
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counter  to  the  direction  of  the  swirl.  The  results  of  these 
tests  provide  design  criteria  for  selecting  injection  nozzle 
configuration  and  estimating  the  required  injection  flow  for 
a turbine  of  given  performance  characteristics.  A potential 
solution  of  the  flow  through  a two-dimensional  radial  cascade 
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A study  was  undertaken  to  evaluate  a proposed  method  of  preventing 
the  occurrence  of  draft  tube  surge  In  hydroelectric  pump-turbines.  It 
Is  generally  accepted  that  the  origin  of  draft  tube  surge  is  related 
to  the  amount  of  angular  momentum  or  swirl  left  in  the  discharge 
flow  from  the  turbine.  The  swirl  gives  rise  to  unstable  flow  patterns 
which  causes  pressure  pulsations  described  as  draft  tube  surge.  The 
adverse  consequences  of  surge  are  noise,  vibration,  vertical  move- 
ment of  the  runner  and  shaft,  variations  in  power  output,  and  pressure 
pulsations  in  the  penstock. 

Realizing  that  the  onset  of  draft  tube  surge  is  associated 
with  the  amount  of  swirl  present  in  the  discharge  flow,  the  obvious 
solution  would  consist  of  eliminating  the  swirl.  Straightening 
vanes  and  fins  located  in  the  draft  tube  have  been  suggested  and 
tried,  but  have  resulted  in  either  efficiency  losses  or  structural 
and  cavitation  damage  as  reported  in  [1].  Injection  of  air  has, 
in  some  cases,  reduced  the  magnitude  of  the  pressure  pulsations. 
Appendages  attached  in  the  draft  tube  such  as  a hollow  cylinder 
to  contain  the  vortex  core  or  solid  fairings  to  reduce  the  intensity 
of  the  vortex  have  met  with  limited  success.  The  majority  of  the 
methods  that  have  been  attempted  cause  either  excessive  energy 
losses,  result  in  cavitation  damage,  or  induce  excessive  structural 
vibrations. 

A method  that  reduces  the  rotation  in  the  draft  tube  flow  at 
off-design  conditions  and  does  not  effect  the  efficiency  of  the 
machine  when  it  is  operating  at  or  near  its  point  of  best 
efficiency  is  necessary.  It  must  not  consist  of  appendages  in 
the  draft  tube  or  impair  the  performance  of  the  machine  when  it 
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is  operating  as  a pump.  The  elimination  of  either  positive  or 
negative  swirl  should  be  provided. 

The  method  that  is  considered  and  evaluated  in  this  report 
consists  of  a series  of  flush  mounted  nozzles  located  in  the  draft 
tube  immediately  downstream  of  the  turbine  rotor.  A schematic  of 
the  arrangement  is  shown  in  Figure  (1).  The  nozzles  will  inject 
fluid  in  the  draft  tube  counter  to  the  peripheral  motion  of  the 
discharge  flow.  Although  not  indicated  in  the  figure,  two  separate 
rows  of  nozzles  could  be  provided  with  one  row  capable  of  injecting 
fluid  in  a direction  opposite  to  that  of  the  other.  An  arrangement 
such  as  this  would  permit  the  reduction  of  either  positive  or  negative 
swirl.  It  is  envisioned  that  the  nozzles  will  be  directly  connected 
to  the  high  pressure  fluid  in  the  spiral  casing  or  penstock.  An 
appropriate  system  of  valving,  activated  on  the  basis  of  wicket 
gate  opening,  would  control  the  flow  to  the  nozzles  and  could  control 
the  number  of  nozzles  discharging. 

An  alternate  arrangement  to  that  shown  in  Figure  (1)  would 
consist  of  swirl  nozzles  individually  connected  to  the  spiral 
casing  by  separate  piping  and  valving.  The  nozzles  would  be  flush 
mounted  on  the  wall  of  the  draft  tube  and  thereby  provide  a minimum 
flow  disturbance  when  the  swirl  nozzles  are  inactive.  Their  presence 
should  not  cause  any  loss  in  efficiency  or  cavitation  resistance 
of  the  machine  when  they  are  inactive. 

Contained  in  this  report  is  an  experimental  study  of  the 
efficiency  of  various  nozzle  geometries  with  respect  to  the  removal 
of  swirl  in  the  draft  tube.  The  resulting  empirical  data  is  applied 
to  specific  turbine  applications  to  illustrate  a method  to  predict 
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the  amount  of  bypass  fluid  and  to  specify  the  nozzle  geometry  required 
to  suppress  draft  tube  surge. 

In  addition  to  the  injection  nozzle  studies  a potential  solution 
of  the  flow  through  a two-dimensional  radial  cascade  of  airfoils  is 
presented,  where  the  airfoils  in  the  cascade  can  be  of  any  arbitrary 
shape.  This  solution  is  applicable  to  the  wicket  gates  of  hydraulic 
turbines  and  was  developed  as  an  aid  in  the  prediction  of  draft  tube 
surge.  The  method  of  solution  consists  of  transforming  the  radial 
cascade  into  a two-dimensional  rectangular  cascade  by  using  a conformal 
transformation.  An  existing  cascade  analysis  program,  known  as  the 
Douglas-Neumann  Cascade  Program,  was  then  employed  to  obtain  the  potential 
flow  solution  in  the  transformed  plane.  Using  a model  of  a wicket 
gate  and  draft  tube  system,  with  air  as  the  fluid,  experiments  measuring 
the  flow  angles  downstream  of  the  wicket  gates  were  conducted  to  evaluate 
the  accuracy  of  the  prediction  method.  Very  good  agreement  was  found 
between  the  measured  and  predicted  fluid  angles. 

2.  INJECTION  NOZZLE  STUDIES 
2.1  Estimate  of  Bypass  Flow 

A theoretical  estimate  of  required  bypass  flow  to  implement  the 
proposed  means  of  surge  elimination  shown  in  Figure  1 shall  be  established 
first.  This  fluid  shall  be  assumed  to  be  piped  from  the  spiral 
casing  and  discharged  into  the  draft  tube.  The  bypassing  of  such  fluid 
represents  an  energy  loss  that,  for  purposes  of  analysis,  will  be 
considered  as  totally  unrecoverable.  On  this  basis,  the  swirl  in  the 
draft  tube  must  be  reduced  to  some  specified  value  with  a minimum  rate 
of  injected  fluid  if  the  described  technique  is  to  be  successful.  An 


-11- 


WSG:AMY:TS: jep 
20  March  1978 


, 


engineering  estimate  of  the  ratio  of  bypassed  flow  to  flow  through 
the  turbine  can  be  obtained  on  the  basis  of  momentum  considerations. 

The  flux  of  angular  momentum  (ft)  and  draft  tube  flow  (Q)  can  be 
determined  at  the  exit  of  the  wicket  gates.  Figure  2,  as: 


and 


Q 


( V cos  a dA  ■ V cos  a 2ir  r B 


ft  - p 


V cos  a (r  V sin  a)  dA 


PQtj 


(1) 

(2) 


The  quantity  B is  wicket  gate  depth  and  Vq  is  the  peripheral  component 
of  velocity  at  wicket  gate  exit. 


} 

. 
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From  equations  (1)  and  (2)  the  momentum  parameter  ^z-  can  be  converted 

V0T  W 

to  a ratio  of  ~ — by  assuming  a free  vortex  flow  in  the  draft  tube,  where 
aT 

the  quantity  Vg^  is  the  peripheral  velocity  at  the  draft  tube  wall,  and 
vAt  is  the  axial  velocity  at  the  entrance  to  the  draft  tube  (or  exit  of  the 
turbine).  Applying  the  law  of  conservation  of  angular  momentum  to  the 
flow  as  it  passes  from  the  wicket  gates  and  enters  the  draft  tube  gives: 


D 

2 


“ ve„ 


%ri 


(3) 


where  D is  the  diameter  of  the  draft  tube  inlet.  By  (1),  (2)  and  (3), 
the  momentum  parameter  becomes: 


S ■ Momentum  Parameter 


Vei  rl  P 2 VqT 

Q " » vAt 


' f f * * 


12- 


WSG:AMY:TS: jep 
20  March  1978 
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(V(T> 


The  ratio  of  peripheral  to  axial  velocity  can  now  be  written  in  terms 
of  the  momentum  parameter: 


V6X 

vAt 


1.57 


1.57  (S) 


(4) 


It  is  proposed  to  reduce  the  swirl  that  exists  in  the  draft  tube  at 
any  given  gate  opening,  Vq^,,  to  that  level  of  swirl  which  exists  prior 
to  the  occurrence  of  periodic  draft  tube  surge,  Vqt»  This  level  of 
swirl  would  be  attained  when  the  angular  momentum  of  the  draft  tube  flow 
and  that  of  the  injected  fluid  is  summed  to  give  a resultant  momentum  para- 
meter, S,  which  is  less  than  that  at  which  periodic  draft  tube  surge  is 

* 

predicted  to  occur,  S . 

ic 

If  the  critical  momentum  parameter  S is  not  to  be  exceeded  in  the 

draft  tube,  the  amount  of  fluid  that  must  be  injected  through  the  swirl 

nozzles  can  be  evaluated.  The  quantities  of  flow  and  velocity  associated 

with  the  critical  momentum  parameter  shall  be  noted  by  an  asterisk,  and 

subscript  (T)  denotes  flow  and  velocities  associated  with  the  draft  tube. 

The  quantity  of  bleed  fluid  required  to  obtain  the  critical  momentum 
* 

parameter  S in  the  draft  tube  can  be  determined  by  equating  the  difference 

between  the  angular  momentum  from  the  turbine  discharge  Q^rVgT  and  that 

of  the  bleed  fluid  Qgr/2gH  to  the  angular  momentum  at  the  critical 
* * 

condition  QTrV0T.  Thus, 


(QTrV0T)  Qfir»^2gH 

* *~7  " , * * 

<QTrVAl)  (QTrVAT) 
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at 


Qj  *^gH 
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Qt  VAn 


- 1.57  S 


The  above  can  be  rearranged  to  provide  a ratio  of  bypass  flow  to  draft 
tube  flow  in  terms  of  the  momentum  parameters  (S)  and  (S  ) : 

* 


1.57  V 


vAt  _ .. 

at  ^ - (vf>  s 1 


at 


at.  *. 
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'2gH 


(5) 


The  decrease  in  efficiency (Ar|)  of  the  turbine  due  to  bleeding  bypass 
fluid  to  the  swirl  nozzles  assuming  that  none  of  the  bypass  fluid  energy 
is  recovered  is: 


\S 


An  - Cioo)  (^) 


(6) 


If  Equation  (5)  is  considered  for  a given  (S),  the  numerator 
is  essentially  constant  and  independent  of  turbine  head.  This  is 
based  on  the  assumption,  for  cavitation  purposes,  that  the  velocity 
at  the  design  condition  in  the  draft  tube  throat  for  conventional 

QB 

turbine  installations  has  some  upper  limit.  The  ratio  of  — specified 

^T 

by  (5)  decreases  as  some  reciprocal  function  as  head  is  increased. 

On  this  basis.  Equation  (6)  indicates  the  decrease  in  efficiency, 
due  to  bleeding  bypass  fluid  to  the  swirl  nozzles,  will  be  less  in 
a high  head  turbine  than  in  a low  head  turbine.  The  performance 
characteristics  for  a turbine  such  as  presented  for  the  Grand  Coulee 
Third  Power  Plant  Units  in  [3]  provides  the  quantities  required  in 
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equation  (S)  to  estimate  the  ratio  of  bypass  to  draft  tube  flow  at  a 
given  momentum  parameter. 

2.2  Test  Apparatus 

Investigations  of  the  performance  of  various  Injection  nozzle  geometries 
were  conducted  In  an  air  flow  facility  which  Is  shown  In  Figure  3. 

Similar  studies  conducted  In  [2]  and  [4]  have  shown  that  the  use  of 
air  as  a working  fluid  has  been  quite  successful  in  characterizing  the 
momentum  parameter  in  draft  tube  flow  and  predicting  the  occurrence 
of  draft  tube  surge.  The  main  air  supply  was  provided  by  a variable 
speed  centrifugal  blower  and  was  measured  by  an  orifice  meter  arrangement. 

The  metered  air  flow  was  diffused  as  it  entered  a stilling  chamber  where 
screens  were  used  to  reduce  flow  turbulence  across  the  stilling  chamber 
flow  field.  The  air  radially  entered  a cylindrical  draft  tube  through 
wicket  gate  type  swirl  vanes.  The  angle  between  the  vanes  and  a radial 
line  could  be  set  at  any  angle  between  0°  (radial  position)  and  82.5° 

(closed  position).  It  can  be  shown  that  the  dimensionless  momentum 
parameter,  , reduces  to  K (tan  a)  for  this  test  facility,  where 
K = constant  and  a is  the  average  flow  angle  leaving  the  wicket  gates. 

Flow  traverse  data  taken  across  the  trailing  edge  of  the  wicket  gates 
using  a prism  probe  were  used  to  determine  flow  angles  for  various  wicket 
gate  settings.  These  data  were  used  to  develop  a relation  for  the  momentum 
parameter  as  a function  of  gate  setting.  Injection  fluid  provided  by 
an  auxiliary  air  supply  and  measured  by  an  orifice  meter  arrangement, 
entered  the  manifold  of  the  injection  nozzle  test  section  located  just 
below  the  wicket  gates.  Any  ratio  of  injection  fluid  flow,  Q , to  fluid 

D 

flow  through  the  gates  (turbine  flow),  Q^,  could  be  provided  by  the 
proper  setting  of  a butterfly  valve  on  the  auxiliary  air  supply  blower. 
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equation  (5)  to  estimate  the  ratio  of  bypass  to  draft  tube  flow  at  a 
given  momentum  parameter. 

2.2  Test  Apparatus 

Investigations  of  the  performance  of  various  Injection  nozzle  geometries 
were  conducted  In  an  air  flow  facility  which  is  shown  in  Figure  3. 

Similar  studies  conducted  in  [2]  and  [4]  have  shown  that  the  use  of 
air  as  a working  fluid  has  been  quite  successful  in  characterizing  the 
momentum  parameter  in  draft  tube  flow  and  predicting  the  occurrence 
of  draft  tube  surge.  The  main  air  supply  was  provided  by  a variable 
speed  centrifugal  blower  and  was  measured  by  an  orifice  meter  arrangement. 

The  metered  air  flow  was  diffused  as  it  entered  a stilling  chamber  where 
screens  were  used  to  reduce  flow  turbulence  across  the  stilling  chamber 
flow  field.  The  air  radially  entered  a cylindrical  draft  tube  through 
wicket  gate  type  swirl  vanes.  The  angle  between  the  vanes  and  a radial 
line  could  be  set  at  any  angle  between  0°  (radial  position)  and  82.5° 

(closed  position).  It  can  be  shown  that  the  dimensionless  momentum 
parameter,  , reduces  to  K (tan  a)  for  this  test  facility,  where 
K - constant  and  o is  the  average  flow  angle  leaving  the  wicket  gates. 

Flow  traverse  data  taken  across  the  trailing  edge  of  the  wicket  gates 
using  a prism  probe  were  used  to  determine  flow  angles  for  various  wicket 
gate  settings.  These  data  were  used  to  develop  a relation  for  the  momentum 
parameter  as  a function  of  gate  setting.  Injection  fluid  provided  by 
an  auxiliary  air  supply  and  measured  by  an  orifice  meter  arrangement, 
entered  the  manifold  of  the  injection  nozzle  test  section  located  just 
below  the  wicket  gates.  Any  ratio  of  Injection  fluid  flow,  Q , to  fluid 
flow  through  the  gates  (turbine  flow),  QT,  could  be  provided  by  the 
proper  setting  of  a butterfly  valve  on  the  auxiliary  air  supply  blower. 
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The  nozzle  test  section  provided  three  points  of  fluid  Injection  equally 
spaced  In  a circumferential  plane  perpendicular  to  the  draft  tube 
centerline.  The  nozzle  geometries  tested  are  shown  In  Table  I.  It 
should  be  pointed  out  that  in  Table  1,  A^,  nozzle  area,  was  defined 
as  the  sum  of  the  three  separate  nozzle  areas.  After  passing  through 
the  nozzle  test  section,  the  air  flow  entered  the  cylindrical  draft 
tube  (L/D~4.56).  Pressure  taps  and  probe  holes  were  located  along  the 
length  of  the  tube.  The  cylindrical  draft  tube  could  be  replaced  with 
an  elbow  type  (Fontenelle)  draft  tube. 

The  unsteady  pressure  produced  by  the  swirling  flow  in  the  draft 
tube  was  monitored  at  the  last  pressure  tap  (L/D=4.40)  on  the  draft 
tube  by  a dynamically  calibrated  differential  pressure  transducer. 

The  pressure  signal  was  analyzed  by  a Spectral-Dynamics  real  time  analyzer 
in  conjunction  with  an  ensemble  averager  from  which  the  frequency,  f, 
of  the  surging  condition  could  be  established. 

Figure  4 shows  typical  data  obtained  from  the  analyzer.  The  root- 
mean-square  value  of  the  unsteady  pressure,  2 , at  the  surging  condition 

was  read  from  an  RMS  meter  after  the  signal  was  passed  through  a line 
filter  with  a characteristic  48  dB  per  octave  roll-off.  The  band  pass 
of  the  filter  was  set  at  20  Hz  to  120  Hz  for  all  tests.  The  frequency 
of  all  observed  pressure  surges  was  between  35  and  70  Hz.  Directly  across 
from  the  pressure  pickup  tap  used  on  the  draft  tube  was  a probe  hole  used 
for  flow  field  measurements.  A prism  probe  was  used  to  obtain  the  time-mean 
total  pressure,  static  pressure  and  fluid  flow  angle  as  a function  of  radial 
distance  across  the  draft  tube  for  various  flow  conditions.  The  basic 
parameters  used  in  this  study  to  describe  the  surging  condition  were  the 
dimensionless  groups  suggested  by  [5];  namely,  the  pressure  parameter 

^ 2 o 

(°  *y2/pQ  ) and  the  frequency  parameter  (fD  /Q). 
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2.3  Discussion  of  Experimental  Results 

The  effect  of  various  fluid  injection  angles  and  nozzle  geometries 
were  investigated  on  the  basis  of  their  effectiveness  in  reducing  the  unsteady 
pressure  amplitude  measured  in  the  draft  tube.  All  tests  were  at  a Reynolds 
Numbers  above  80,000  as  recommended  in  [4]  and  at  velocities  low  enough 
to  prevent  compressible  effects.  A typical  spectral  analysis  Indicating 
the  reduction  in  the  amplitude  of  the  unsteady  surge  pressure  is  shown 
by  Figure  (4).  It  is  apparent  that  both  the  amplitude  and  the  frequency 
of  the  surge  are  reduced  as  the  ratio  of  Q„/Q  is  increased.  The  spectral 
analysis  indicates  no  pressure  peaks  at  a Q_/Q_  of  16%  indicating  that 
surge  has  been  eliminated  at  this  condition  of  fluid  injection. 

The  angle  of  fluid  injection,  as  defined  in  Table  I,  was  studied  with 
respect  to  its  effect  on  reducing  surge  for  a given  swirl  in  the  draft 
tube.  Four  angles  of  injection  were  investigated,  these  being  0°,  15°, 

30°  and  45°  as  defined  by  Table  I.  The  nozzle  geometry  used  in  these 
tests  was  # 5,  described  in  Table  I,  and  the  results  are  shown  in  Figure 
(5)  for  a momentum  parameter  of  0.8  and  1.18.  The  most  significant  result 
being  that,  up  to  the  45°  angle  tested,  the  reduction  of  surge  pressure 
amplitude  is  independent  of  injection  angle.  The  lower  end  points  of 
the  curves  in  Figure  (5)  represent  that  ratio  of  Q_/Q_  required  to  completely 
eliminate  surge  and  it  is  evident  that  the  quantity  of  Qg/Qj  required, 
is  equal  for  the  four  injection  angles  tested. 

An  Injection  angle  of  30°  was  selected  for  the  tests  to  investigate 
the  effect  of  the  length  to  width  ratio  of  the  Injection  nozzle  as  well 
as  the  ratio  of  the  area  of  the  injection  nozzles  to  the  draft  tube  throat 
area.  The  ratio  Qg/Qj  is  plotted  against  the  surge  pressure  parameter 
for  a series  of  nozzle  geometries  for  momentum  parameters  of  0.8  and  1.18 
in  Figure  6.  The  lower  extremity  of  each  curve  represents  a condition 
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where  surge  Is  eliminated  and  thereby  indicates  the  ratio  of  QD/Q_ 

required  to  eliminate  surge  for  a given  A^/AT  ratio.  Figure  (6)  illustrates 

that  for  a given  ratio  of  Qg/Q^t  the  reduction  in  the  surge  pressure 

amplitude  is  greater  as  A^/^  decreases.  This  results  directly  from 

the  fact  that  to  inject  a given  ratio  of  Q /Q  into  the  draft  tube,  the 

15  1 

head  and,  hence,  the  velocity  must  increase  as  A^/A^,  decreases.  Since 
the  momentum  of  the  injected  fluid  increases  as  the  square  of  the  spouting 
velocity  from  the  nozzles,  smaller  ratios  of  Q /Q  are  required  to 
eliminate  surge  as  the  head  on  the  turbine  increases. 

Table  I lists  the  various  nozzle  geometries  tested  and  defines 
their  characteristic  dimensions.  The  rather  surprising  result  being 
that  for  the  various  L/W  nozzles  tested,  no  measurable  difference  was 
observed  in  the  ratio  of  Q_/Q_  required  to  eliminate  surge  at  a given 
Ajj/A^,.  This  result  suggests  that  the  choice  of  a nozzle  geometry  to 
obtain  a given  A^/A^,  would  consist  of  selecting  that  geometry  which 
would  give  the  least  hydraulic  losses  and  permit  the  easiest  fabrication. 

On  this  basis,  a cylindrical  nozzle  which  would  be  elliptical  at  its 
intersection  with  the  barrel  of  the  draft  tube  would  seem  appropriate. 

The  tangential  and  axial  velocity  components  of  the  flow  in  the 
draft  tube  were  measured  at  a station  of  L/D=4.40,  indicated  on 
Figure  (3),  by  means  of  a prism  probe  which  provided  a time-i-veraged 
reading  of  the  local  static,  total  pressure  and  flow  angularity. 

The  velocity  distributions  are  shown  by  Figure  (7)  for  momentum 
parameters  of  0.24,  0.41,  0.80  and  1.18.  The  dashed  1 >e  in  these 
figures  represents  similar  traverse  data  having  the  <•>  cket  gates 
at  a setting  corresponding  to  a momentum  parameter  o.  .8  but  with 
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fluid  injected  to  eliminate  surge.  It  is  evident  that  the  reversal 
of  axial  flow  at  the  axis  of  rotation  is  eliminated  and  considerable 
swirl  has  been  removed  from  the  flow  when  fluid  is  injected. 

Without  injection,  the  outer  portion  of  the  flow  approaches 
a free  vortex,  or  constant  angular  momentum,  distribution  with 
solid  body  rotation  near  the  center.  The  peripheral  velocity 
distribution  downstream  of  a turbine  runner  may  deviate  from  that 
indicated  in  Figure  (7).  The  effect  of  such  variation  has  not 
been  investigated  in  the  present  studies,  but  deserves  future 
consideration. 

2.4  Comparison  of  Test  Results  with  Predictions 

The  ratio  of  bleed  flow  to  turbine  flow  was  previously  derived 
in  Equation  (5),  assuming  ideal  momentum  transfer  between  the 
injected  and  draft  tube  flow.  On  the  basis  of  the  preceding 
experimental  data,  it  is  possible  to  obtain  an  empirical  relation 
that  applies  to  pump-turbines.  The  bypass  flow  can  be  expressed 
as 


Qb  = K1(2g  H)1/2  Ajj  , 

where  is  a nozzle  coefficient  and  assumed  equal  to  0.9.  The 
flow  through  the  turbine  can  be  written  as 


(7) 


(8) 
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The  ratio  of  bleed  flow  to  turbine  flow  is  then 


Qb/Qt  - 39,100  n H3/2 


|(451  n^^H3/2) 


• (9) 


English  Units 

The  results  of  experimental  data  indicating,  for  a given  A^/A^,,  the 
ratio  of  Qg/Q^,  required  to  eliminate  surge  at  two  momentum  parameters 
are  shown  in  Figure  (8).  A curve  fitting  of  this  data  indicates 
a relation 


Qb/Qt  - 0.719  (An/At) 


0.401 


(10) 


for  S=1.18  and. 


Qb/Qt  = 0.473  (An/At) 


0.366 


Ul) 


for  S=0.80. 

Considering  a momentum  parameter  of  1.18,  it  is  possible  by 


substituting  (10)  in  (9)  to  obtain  a relation  for  Qg/QT  of, 


104.7  t(n)  H 

^ ' English  Units 

In  similar  fashion,  the  substitution  of  (11)  in  (9)  provides  a 


relation  of  QR/QT  for  a momentum  parameter  of  0.8 
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qb/qt 


1460  (V1)0'578  H°'867 


111.5  (1^,0-578  „0.867 


. (13) 


v -'English  Units 

An  elbow-type  draft  tube  was  also  evaluated  to  Indicate  the 

Influence  of  draft  tube  shape.  The  draft  tube  tested  was  similar 

to  the  Fontonelle  configuration  described  in  [4] . The  data  and 

equations  relating  A1J/AT  to  Qfi/QT  required  to  eliminate  surge  for 

this  draft  tube  are  shown  in  Figure  9.  By  applying  the  same  analysis 

as  indicated  for  the  cylindrical  draft  tube,  the  ratio  of  bypass 

flow  to  turbine  flow  for  the  elbow-type  draft  tube  was  calculated. 

Using  performance,  geometrical  data  and  a speed  coefficient  of 

<J>=0.8  for  the  model  of  the  Grand  Coulee  Third  Power  Plant  Units 

given  in  [3],  the  bypass  flow  rates  were  determined  for  a cylindrical, 

and  an  elbow-type  draft  tube.  The  results  are  shown  in  Table  II. 

Table  II 

GRAND  COULEE  THIRD  POWER  PLANT  BYPASS  RATES 


S=0.8 


S-1.18 


Predicted 
Ideal 
(Eq.  (5)) 

Predicted 

Empirically 


Qb/Qt=4.49% 

Qb/Qt=10.3% 

Qb/Qt=11.5% 


Qb/Qt=3.97% 

Qb/Qt*13.9%  (cylindrical) 
Qb/Qt=13.0%  (elbow) 


The  comparison  indicates  that  the  empirically  derived  quantity  of 
Qb/Qt  required  to  eliminate  surge  at  a momentum  parameter  of  1.18 
is  about  3.5  times  greater  than  that  predicted  assuming  an  ideal 
momentum  transfer.  At  the  lower  momentum  parameter,  the  difference 
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between  the  estimated  values  of  Qg/QT  Is  not  as  great.  The  difference 
between  predicted  and  measured  values  of  Qg/Q^<  required  to  prevent  draft 
tube  surge  Is  preimary  due  to  the  initial  assumptions  that  the  peripheral 
velocity  distribution  in  the  draft  tube  was  a potential  vortex  and  that  the 
axial  velocity  was  uniform.  It  is  apparent  from  Figure  7 that  these 
assumptions  are  not  realistic  when  surge  occurs  in  the  draft  tube. 

It  should  be  emphasized  that  solution  of  Equation  (12)  and 
(13)  in  conjunction  with  Figure  (8)  also  specifies  the  A^/A^, 
required,  for  a given  turbine,  to  eliminate  surge  at  a value  of  S 
equal  to  0.8  or  1.18. 

2.5  Summary  and  Recommendations 


Experimental  studies  with  respect  to  eliminating  draft  tube  surge 
by  injecting  fluid  into  the  draft  tube  counter  to  the  existing  swirl 
resulted  in  the  following  observations. 

As  the  ratio  of  injected  fluid  Q /Q  is  increased,  both  the  surge 
frequency  and  the  unsteady  surge  pressure  decreases  in  magnitude. 

The  use  of  various  nozzles  of  a fixed  A^A^,  indicates  no  measurable 
change  in  the  effectiveness  of  the  nozzle  in  reducing  or  eliminating 
surge.  Similarly,  nozzle  injection  angles  up  to  45°  show  no  difference 
in  their  effectiveness  in  reducing  surge.  On  this  basis,  the  nozzle 
geometry  would  be  selected  primarily  on  the  basis  of  hydraulic  efficiency 
and  ease  of  fabrication. 

Empirical  data  are  presented  which  permits  estimating  the  area 
of  the  injection  nozzles  and  the  quantity  of  bleed  fluid  required 
for  turbines  of  specified  performance  characteristics. 

At  a given  momentum  parameter,  the  quantity  Q /Q  required  to 

B T 

eliminate  surge  decreases  as  higher  head  turbines  or  lower  specific 
speed  machines  are  considered. 
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Further  Investigations  are  required  to  gather  data  indicating 
the  influence  of  draft  tube  shape  on  the  effectiveness  of  fluid 
injection  in  reducing  draft  tube  surge. 

The  hydrodynamic  effects  such  as  cavitation  arising  from 
the  interaction  of  the  draft  tube  flow  with  the  high  velocity 
injection  jets  should  be  experimentally  investigated. 

3.  ANALYSIS  OF  FLOW  FIELD  THROUGH  WICKET  GATE 
3.1  Preliminary  Assumptions  and  Governing  Equations 

Experimental  studies  have  shown  that,  for  a given  draft  tube 

shape,  the  frequency  and  pressure  parameters  are  independent  of 

Reynolds  number  for  Reynolds  numbers  greater  than  1 x 10"*  [2,5,6]. 

This  is  an  important  consideration,  since  turbine  prototype  Reynolds 

numbers  exceed  1 x 10^,  and,  thus,  studies  of  surge  can  be  conducted 

eliminating  Re  as  a variable.  It  has  also  been  found  experimentally 

[2,5,6]  that,  for  a particular  draft  tube,  the  frequency  and  pressure 

parameters  correlate  with  the  momentum  parameter,  and  a critical 

2 

value  of  the  momentum  parameter  (QD/pQ  ) exists  above  which  surge 
will  occur.  These  facts  have  enabled  extensive  studies  to  be  conducted 
on  models  of  various  types  of  draft  tubes  and  for  different  values 
of  (L/D) , from  which  the  pressure  and  frequency  parameters  can  be 
correlated  with  the  momentum  parameter  and  used  to  predict  surge  for 
prototype  installations  [A]. 

The  correlation  of  draft  tube  surge  characteristics  with 
the  momentum  parameter  brings  the  realization  of  the  importance 
of  accurately  calculating  this  parameter.  The  most  convenient 
way  to  evaluate  the  momentum  parameter  is  to  determine  the 
angular  momentum  of  the  flow  leaving  the  wicket  gates  and  subtract 
from  this  value  the  angular  momentum  removed  by  the  turbine 
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runner.  The  momentum  parameter  for  the  flow  in  the  draft  tube 
can  thus  be  written  as: 


where 


_ PD 
pwQ2 

i 

P * turbine  power  output 
u * angular  velocity  of  the  turbine. 


flD 

PQ2 


QD 


PQ 


(14) 


The  second  term  on  the  right  of  Equation  (14)  represents  the 
change  in  angular  momentum  across  the  turbine  runner  and  is 
evaluated  from  the  turbine  performance  characteristics  usually 
obtained  from  model  studies.  The  first  term  on  the  right  is  the 
momentum  parameter  for  the  flow  leaving  the  wicket  gates  and  is 
currently  evaluated  using  the  graphical  approach  of  [5]  and  illustrated 
by  Figure  10.  With  this  approach,  the  flow  is  assumed  to  leave 
the  wicket  gates  perpendicular  to  the  minimum  cross  section 
between  the  trailing  edge  of  one  wicket  gate  and  the  adjacent 
wicket  gate.  Referring  to  the  nomenclature  defined  in  Figure  10, 
the  momentum  parameter  at  the  exit  of  the  wicket  gates  can  be  computed 
by  the  following  expression: 


ftD 


PQ 


D r^  sin  a 
BNs 


(15) 


where 


B * depth  of  the  wicket  gates 
N = number  of  wicket  gates. 
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Although  the  most  convenient  way  of  calculating  the  momentum 
parameter  Is  to  use  the  approach  expressed  by  Equation  (14),  it 
is  mentioned  in  References  [3]  and  [4]  that  the  graphical 
method  of  evaluating  the  first  term  of  Equation  (14)  may  be 
seriously  affecting  the  accuracy  of  using  the  momentum  parameter 
to  correlate  experimental  data  and  to  predict  the  occurrence 
and  characteristics  of  draft  tube  surge.  For  this  reason,  it  is 
felt  that  more  attention  must  be  given  to  the  details  of  the  flow 
through  the  wicket  gates.  It  is  also  realized  that  the  upstream 
stay  vanes  or  the  inlet  spiral  may  influence  the  flow  leaving 
the  wicket  gates;  however,  these  factors  obviously  cannot  be 
accounted  for  until  an  analysis  is  employed  which  considers  the 
details  of  the  actual  flow  process,  which  is  beyond  the  capability 
of  the  graphical  approach. 

With  the  above  considerations  in  mind,  the  following 
presents  a method  to  obtain  a potential  solution  of  the  flow 
through  the  wicket  gates  and  presents  the  results  of  an 
experimental  study  conducted  to  evaluate  the  accuracy  of  this 
approach.  A potential  flow  solution  was  developed,  because  it 
represents  the  first  step  usually  taken  in  a problem  of  this 
nature,  which  can  be  used  as  a basis  for  a more  exact  model,  if 
necessary.  Since  the  flow  through  the  wicket  gates  is  an  accelerating 
flow,  the  potential  flow  solution  is  expected  to  yield  quite 
satisfactory  results  by  itself. 

Several  methods  for  analyzing  the  potential  flow  through 
two-dimensional  rectangular  cascades  are  presently  available. 

One  of  the  most  general  of  these,  capable  of  handling 
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airfoils  of  any  arbitrary  shape.  Is  a method  developed  at 
Douglas  Aircraft  Corporation  and  referred  to  as  the  Douglas- 
Neumann  Cascade  program.  A conformal  transformation  of  a radial 
flow  cascade  to  a rectangular  cascade  enables  the  use  of  the 
Douglas-Neumann  Cascade  program  to  analyze  the  flow  in  the 
transformed  plane. 

The  first  assumption  which  is  necessary  to  analyze  the 
flow  through  the  wicket  gates  is  that  the  flow  is  two-dimensional. 

For  most  turbine  installations  this  is  a very  reasonable  assump- 
tion, since  the  wicket  gates  themselves  are  two-dimensional,  and 
the  flow  at  the  inlet  and  exit  to  the  wicket  gates  should  be 
predominately  two-dimensional.  A schematic  of  a turbine  cross 
section  is  presented  in  Figure  11,  showing  the  inlet  spiral,  stay 
vanes,  wicket  gates,  turbine  runner  and  draft  tube.  This  figure 
is  provided  to  illustrate  that  the  flow  passage  in  the  region  of 
the  wicket  gates  is  usually  a straight  section,  which  is  the 
justification  for  assuming  the  flow  is  two-dimensional.  With  the 
two-dimensional  flow  assumption,  the  geometry  of  a segment  of  the 
wicket  gate  system  to  be  analyzed  is  presented  in  Figure  12.  This 
, figure  also  indicates  the  coordinate  system  which  is  used. 

The  other  necessary  assumptions  are  those  which  are  required 
to  enable  the  flow  to  be  considered  a potential  flow.  It  is, 
therefore,  assumed  there  are  no  body  forces  and  that  the  flow  is 
steady,  incompressible,  inviscid  and  irrotational.  The  lnviscid 
flow  assumption  is  the  only  assumption  requiring  some  justifica- 
tion. Because  the  flow  is  accelerating  through  the  wicket 
gates,  the  effects  of  the  boundary  layers  will  be  minimal,  and 
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thus,  an  lnvlscid  solution  is  actually  quite  realistic  in  this 


situation 


With  the  above  assumptions,  the  governing  equations  for 


the  flow  are  the  equation  for  zero  vorticity  and  the  continuity 


equation.  These  equations  expressed  in  polar  coordinates 


and  in  terms  of  dimensionless  variables  are  as  follows 


continuity 


A stream  function  (i|>)  is  defined,  such  that  by  its  definition 


it  satisfies  the  continuity  equation 


Substituting  Equations  (18)  and  (19)  into  Equation  (16)  yields  the 


familiar  Laplace  equation  in  polar  coordinates 


Equation  (20)  is  solved  indirectly  through  a transformation  to 


obtain  the  flow  field  solution  through  the  wicket  gates.  The 


boundary  conditions  which  must  be  satisfied  in  conjunction  with 


Equation  (20)  are 
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1.  A specified  flow  direction  at  infinity  or  far  upstream. 

2.  The  wicket  gates  must  constitute  a streamline,  or  in 
other  words,  the  velocity  normal  to  the  wicket  gate 
surface  must  be  zero. 


3.  The  Kutta  condition  must  be  satisfied  at  the  trailing 
edge  of  the  wicket  gates. 

3.2  Transformation  of  the  Flow  Field 

A conformal  transformation  was  used  to  transform  the  wicket 
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gate  system  into  a two-dimensional  linear  cascade,  so  that  the 
flow  in  the  transformed  plane  could  be  analyzed  using  available 
procedures.  Since  the  flow  in  the  original  plane  is  considered 
a potential  flow  and  the  transformation  is  conformal,  the  flow 
in  the  transformed  plane  is  also  a potential  flow  and  can  be 
analyzed  accordingly.  The  transformation  used  is  as  follows: 

Y = 0 (21) 

X *=  In  R . (22) 

By  applying  the  chain  rule  and  Equations  (21)  and  (22),  it 
is  found  that  Equation  (20)  is  transformed  into  Laplace's 
equation  in  cartesian  coordinates.  Equation  (23). 


where 
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velocity  component  In  the  X direction 


(24) 


velocity  component  in  the  Y direction  - 


(25) 


The  transformation  of  Equation  (20)  into  Equation  (23)  was 
carried  out  to  illustrate  that,  with  the  transformation  employed, 
the  flow  in  the  transformed  plane  is  indeed  a potential  flow  and, 
thus,  obeys  the  same  fundamental  laws  as  the  flow  in  the 
original  plane.  In  the  actual  flow  analysis,  the  wicket  gates 
were  transformed,  and  the  linear  cascade  of  airfoils  obtained 
were  analyzed  as  if  this  plane  was  the  real  plane.  The  trans- 
formation is  illustrated  by  Figure  13,  where  the  wicket  gate 
geometry  previously  shown  in  Figure  12  is  shown  in  the  transformed 
plane.  The  circular  arcs  labeled  1,  2,  3 and  4 in  Figure  12 
become  vertical  lines  of  constant  X value  in  Figure  13  and  are 
similarly  labeled  for  a comparison  of  the  two  planes.  It  is 
interesting  to  note  that  the  arc  of  R=1  becomes  the  Y-axis  of 
the  cartesian  coordinate  system  in  the  transformed  plane,  while 
the  point  at  R*0  corresponds  to  X=-°°  in  the  new  plane.  Similarly, 
radial  lines  of  the  original  plane  become  lines  of  constant  Y 
value  in  the  transformed  plane.  The  X-axis  corresponds  to  the 
radial  line  at  6*0. 

The  boundary  conditions  in  the  transformed  plane  are  identi- 
cal to  those  in  the  real  plane.  Since  flow  angles  are  not  changed 
by  a conformal  transformation,  the  specified  inlet  flow  angle 
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is  the  same  in  both  planes.  In  addition,  the  velocity  component 
normal  to  the  surface  of  the  transformed  airfoil  must  be  zero, 
and  the  Kutta  condition  must  be  satisfied,  as  in  the  original 
plane. 

Relationships  between  the  velocity  components  in  the  two 
planes  can  be  found  by  applying  the  chain  rule  and  Equations  (21) 
and  (22)  to  the  derivatives  of  ij>  with  respect  to  0 and  R.  It 
is  first  found  that: 


30  3Y 


(26) 


and 


I li 

3R  R 3X 


(27) 


From  Equations  (26)  and  (27)  and  the  definitions  of  the  stream 
function,  the  desired  relations  between  the  velocities  in  the 
two  planes  are  obtained: 


V * — V 
R R X 


(28) 


V9  - i V,  . (29) 

The  relationship  between  the  magnitudes  of  the  total  velocity 
vectors  can  then  be  expressed  as: 
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original  plane 


Vvn2+ve2  ' i~\F: 


2 2 
+V 

X Y 


(30) 


— V 

R transformed  plane 


3.3  Method  of  Solution  In  the  Transformed  Plane 

The  Douglas  Neumann  cascade  solution  [7]  was  used  for  analyzing 
the  flow  through  the  transformed  cascade  because  this  method 
of  solution  is  capable  of  handling  infinite  cascades  of  airfoils 
of  an  arbitrary  shape.  A solution  is  obtained  by  applying  a 
distribution  of  sources  on  the  surface  of  the  airfoils  in  the 
cascade,  such  that  the  combination  of  the  onset  velocity,  the 
source  at  the  particular  point,  and  the  induced  velocity  from 
the  remaining  source  distribution  satisfy  the  boundary  condition 
of  zero  normal  velocity  at  the  surface.  In  a similar  manner, 
the  vorticity  distribution  on  the  airfoils  is  obtained  by 
employing  the  same  set  of  equations  with  the  velocity  vector  of 
the  sources  rotated  90°. 

The  source  distribution  is  determined  by  representing  the 
airfoils  as  a series  of  straight  line  segments,  with  the  source 
strength  assumed  constant  over  each  segment  and  the  boundary  condition 
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satisfied  at  the  segment  midpoint.  With  this  approximation,  the 
induced  velocity  at  the  midpoint  of  a particular  segment  can  be 
represented  by  a summation  of  a set  of  integrals  representing  the 
induced  velocity  from  the  remaining  segments.  Since  the  source 
strength  is  assumed  constant  over  the  segment,  the  integrals  can 
be  evaluated  analytically,  which  results  in  a set  of  algebraic 
equations  which  must  be  solved  simultaneously  for  the  source 
strength  of  each  segment.  This  technique  will  approach  an  exact 
solution  as  the  number  of  segments  approaches  infinity,  but 
sufficient  accuracy  is  achieved  as  long  as  the  straight  line 
segments  are  small  enough  to  adequately  describe  the  airfoil 
shape. 

The  solution  to  the  general  problem  is  obtained  by 
calculating  the  potential  flow  for  three  basic  flows  consisting 
of  the  solution  for  a flow  with  zero  angle  of  attack,  90°  angle 
of  attack,  and  a pure  circulatory  flow.  These  three  solutions, 
which  all  satisfy  the  boundary  condition  of  a zero  velocity 
component  normal  to  the  surface,  are  then  combined  in  such  a 
manner  to  satisfy  the  specified  inlet  angle  and  the  Kutta  condi- 
tion. The  complete  solution  enables  the  velocity  components 
and  the  static  pressure  to  be  determined  at  any  point  in  the  flow 
field  and  also  yields  the  overall  angle  by  which  the  cascade  turns 
the  flow. 

3. A Discussion  of  Normalizing  Parameters  and  the  Reverse  Transformation 
Thus  far,  the  equations  governing  the  flow  in  both  planes  have  been 
presented,  and  the  method  of  obtaining  a solution  in  the  transformed 
plane  has  been  discussed.  The  equations  relating  the  velocity 
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components  in  the  two  planes  have  also  been  given,  however,  before 

these  equations  can  be  specifically  applied,  consistent  normalizing 

velocities  must  be  selected  in  each  plane.  Since  It  is  known  from  the 

conservation  of  mass  equation  that  the  average  x component  of  velocity 

is  constant  both  upstream  and  downstream  of  a rectangular  cascade,  v^ 

represents  a convenient  parameter  to  normalize  the  velocity  in  the 

transformed  plane.  Although  the  Douglas-Neumann  cascade  solution  uses 

as  a normalizing  parameter  the  modulus  of  the  average  velocity  vector 

upstream  and  downstream  of  the  cascade,  the  output  can  easily  be 

converted  so  that  the  velocity  is  normalized  by  v^. 

In  the  original  plane  the  velocity  component  corresponding  to 

vx  is  the  radial  component.  The  average  radial  velocity  component 

is  not  constant,  but,  once  again  from  the  conservation  of  mass  equation, 

it  is  known  that  r v^  = constant.  Using  this  relationship,  if  a 

reference  value  of  r is  selected,  a corresponding  reference  value  of 

v can  be  obtained.  The  value  of  r which  is  convenient  for  this 
r 

purpose  is  the  radius  of  the  wicket  gate  spindle  centers  (r  ) . Selecting 

s 

rg  as  the  reference  radius  results  in  the  following  definition  of  the 
normalizing  velocity  in  the  original  plane. 

vr  = Q/(2n  B r ) (31) 

s s 

It  is  apparent  from  the  definition  that  vr  represents  the  average 

s 

radial  velocity  which  would  exist  at  the  spindle  radius  if  the 
wicket  gates  were  not  present. 

Rewriting  Equation  (28)  in  terms  of  dimensional  quantities 
allows  the  relationship  between  the  velocities  in  the  planes  to  be 
determined  when  they  are  normalized  using  the  above  parameters.  For 
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the  present,  will  be  normalized  by  the  temporary  variable,  v^,  ^ 


v 

c x 


Vr  r V 


To  relate  the  two  normalizing  velocities,  it  is  known  that  vx  should 

equal  v when  v = vr  and  r-r  . Solving  Equation  (32)  for  v - after 
x * s s ret 

inserting  these  values  yields 


v , = v — 
ref  x r 


Substituting  Equation  (33)  back  into  Equation  (32)  provides  the  new 
relationship  for  relating  the  dimensionless  radial  velocity  in  the 
original  plane  to  the  x component  of  velocity  in  the  transformed  plane. 


Similarly, 


r v 
s x 


R v 

S X 


r — c 
v — 
x r 


r - R - 
v v 

X X 


vQ  R v 

Va  « — 2.  _JL 

0 vr  R — 

s v 

x 


original  plane 


Rs  transformed  plane 


These  equations  provide  the  necessary  relationships  for  the 
reverse  transformation  of  the  solution  from  the  transformed  plane  to 
the  original  plane.  It  should  be  noted  that  the  fluid  angles  are  not 
changed  by  the  transformation,  thus,  the  above  equations  need  only 
to  be  applied  when  interested  in  the  velocity  magnitudes.  Equation 
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(36)  Is  used  in  the  computer  program  to  transform  the  local  velocity 
on  the  surface  of  the  blades,  thus,  enabling  the  pressure  distribution 
on  the  wicket  gate  to  be  determine. 

3.5  Analytical  Study  Conducted 

The  method  of  transforming  the  wicket  gates  and  employing 
the  rectangular  cascade  analysis  to  obtain  a potential  flow 
solution  was  used  to  analyze  the  flow  through  the  wicket  gates 
of  an  air  model.  This  model,  which  was  primarily  fabricated 
to  conduct  draft  tube  surge  studies,  will  be  described  in  the 
next  section  on  the  experimental  study.  The  wicket  gates  shown 
in  Figure  12  are  a scale  drawing  of  the  real  wicket  gates  in  the 
model.  These  wicket  gates  have  no  camber,  a chord  length  of 
3.97  cm  (1.562  in.),  and  a maximum  thickness  of  0.544  cm  (0.214  in.). 

Solutions  were  obtained  for  the  wicket  gates  with  angular 
settings  between  0 and  80°  in  10°  increments.  These  angles  are 
measured  between  the  wicket  gate  chord  and  a radial  line.  The 
closed  position  is  at  the  angular  setting  of  83°.  Two  solutions 
were  obtained  for  each  wicket  gate  setting,  one  with  the  flow 
entering  radially  and  the  other  solution  with  the  flow  at  a 
zero  angle  of  attack  in  the  transformed  plane.  These  two 
solutions  were  obtained  to  evaluate  the  influence  of  the  upstream 
flow  on  the  fluid  exit  angle.  It  was  found  that  for  the  present 
wicket  gate  geometry  the  inlet  angle  did  not  affect  the  exit 
angle.  This  conclusion  cannot  be  made  for  all  systems,  however, 
since  even  for  a potential  flow  solution  the  spacing  between  gates 
will  affect  the  amount  of  turning  the  cascade  can  perform.  The 


-36- 


WSG:AMY:TS:jep 
20  March  1978 


; 


* 

present  method  of  solution  could  be  used  by  testing  various  spaclngs 
to  determine  what  spacing  Is  needed  to  eliminate  upstream  effects. 

The  graphical  method  Indicated  by  Figure  10  was  also  carried 
out,  and  these  results  will  be  presented  along  with  the  experimental 
results  and  potential  flow  solution.  Initial  results  comparing  the 
fluid  exit  angles  predicted  by  the  graphical  approach  and  the  potential 
flow  solution  revealed  very  little  difference.  The  question  was  then 
raised,  would  cambered  wicket  gates  demonstrate  a larger  deviation  between 
the  fluid  angles  predicted  using  the  two  methods?  To  study  this 
possibility,  an  analytical  study  was  conducted  for  cambered  wicket 
gates  similar  to  the  study  for  the  symmetrical  wicket  gates.  For  this 
case,  the  results  from  the  potential  solution  also  revealed  no  significant 
dependence  on  the  inlet  angle.  However,  a slightly  larger  deviation 
was  found  between  the  fluid  exit  angle  predicted  by  the  potential 
solution  and  the  graphical  approach.  The  actual  numerical  results 
will  be  presented  in  a following  section  on  analytical  and  experimental 
results. 

3.6  Experimental  Program 

A sketch  of  the  test  facility  used  to  obtain  experimental 
data  for  comparison  with  the  potential  flow  solution  is  shown  in 
Figure  14.  In  this  facility,  flow  surveys  were  made  at  a 
constant  radius  behind  several  wicket  gate  channels,  with 
measurements  made  every  1.3  degrees.  In  addition  to  the  open 
inlet  configuration  shown  in  the  figure,  tests  were  also  conducted 
with  an  inlet  spiral  installed.  The  tests  with  and  without  the 
inlet  spiral  enabled  the  effects  of  the  inlet  flow  angle  on  the 
exit  flow  from  the  wicket  gates  to  be  investigated.  The  draft 
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tube  has  an  Inlet  diameter  of  15.56  cm  (6.125  In.)  and  the  probe 
was  located  at  a diameter  of  22.86  cm  (9.0  In.). 

Total  pressure,  static  pressure  and  the  flow  angles  were 
measured  using  a three-hole  prism  type  probe.  The  flow  angle 
was  determined  by  nulling  the  probe.  The  total  pressure 
was  measured  from  the  center  hole,  and  the  static  pressure  was 
calculated  from  the  average  pressure  of  the  two  side  holes  and 
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Equation  (37)  was  used  for  calculating  a from  the  momentum 
parameter  obtained  experimentally.  An  average  flow  angle  was  also 
calculated  by  simply  computing  the  arithmetic  average  of  all  the 
measured  local  flow  angles.  These  two  average  angles  were  almost 
identical,  however,  and  future  references  to  the  average  angle  will 
not  mention  the  method  of  calculation. 

3.7  Discussion  of  Analytical  and  Experimental  Results 

The  downstream  flow  angle  obtained  from  the  potential  flow 
solution  in  the  transformed  plane  is  the  angle  of  the  uniform 
flow  theoretically  achieved  at  negative  infinity.  For  all 
practical  purposes,  however,  uniform  flow  of  constant  angularity 
is  achieved  a short  distance  from  the  trailing  edge  of  the 
blades.  When  applying  the  potential  flow  solution  to  a real  situation, 
the  flow  may  or  may  not  become  uniform  before  it  reaches  the  location 
of  the  turbine  runner.  It  should  be  pointed  out  that  even  if  sufficient 
spacing  for  the  flow  to  become  uniform  is  not  available,  the  uniform 
flow  angle  still  applies  for  calculating  the  momentum  parameter, 
since  no  change  in  the  momentum  parameter  will  occur  between  the 
trailing  edge  and  the  location  where  uniform  flow  would  be  achieved. 

With  this  in  mind,  the  flow  angle  calculated  by  the  potential  theory 
for  the  symmetrical  wicket  gate  is  presented  in  Figure  15  as  a function 
of  wicket  gate  angle.  The  angles  from  the  potential  theory  are 
represented  by  the  solid  line. 

The  average  fluid  angles  obtained  experimentally  are  also 
presented  in  Figure  15.  No  difference  was  found  in  the  measured 
angles  for  the  two  sets  of  tests  conducted  both  with  and  without 
the  inlet  spiral,  thus,  only  one  set  of  data  is  presented.  The 
fact  that  the  exit  flow  angles  did  not  change  with  the  different 
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inlet  configurations  is  consistent  with  the  predicted  data  which 
also  showed  no  effect  with  the  change  in  the  angle  of  attack. 

The  graphical  method  for  predicting  the  momentum  parameter. 
Indicated  by  Figure  10,  was  used  to  calculate  the  exit  flow  angle 
from  the  wicket  gates  to  provide  an  additional  comparison  and 
demonstrate  if  any  improvement  is  achieved  by  the  potential 
flow  solution.  The  flow  angle  (a)  defined  in  Figure  10  is  not 
the  average  flow  angle,  however,  since  the  angle  is  defined  within 
the  wicket  gate  passage  and  does  not  represent  the  angle  obtained 
when  the  flow  becomes  uniform.  The  angle  (a)  cai  be  used  with 
Equation  (IS)  to  calculate  the  momentum  parameter,  and  then  with 
the  momentum  parameter  the  average  flow  angle  (a)  can  be  calculated 
from  Equation  (37).  The  average  flow  angles  calculated  in  this 
manner  are  the  data  represented  by  the  circles  in  Figure  15. 

For  the  system  with  the  symmetrical  wicket  gates,  as  shown  by 
Figure  15,  the  results  from  both  the  potential  flow  solution  and 
graphical  method  agree  very  closely  with  the  experimental  data. 

At -the  higher  wicket  gate  angles  the  potential  solution  is  only 
approximately  1 degree  low  in  its  prediction  of  the  fluid  angle. 

At  the  lower  wicket  gate  settings  the  potential  solution  predicts 
values  approximately  2 degrees  higher  than  the  measured  angles. 

On  the  other  hand,  the  graphical  method  shows  its  greatest  deviation 
at  the  higher  wicket  gate  angles,  where  flow  angle  is  under-predicted 
by  approximately  4 degrees.  At  the  lower  wicket  gate  angles,  the 
fluid  angles  from  the  graphical  method  and  the  experimental  data  are 
approximately  equal. 
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With  the  deviation  between  the  predicted  and  experimental 
data  of  such  a small  magnitude,  it  is  difficult  to  conclude 
whether  the  deviation  is  due  to  a minor  shortcoming  of  the 
prediction  methods  or  due  to  an  error  in  the  experimental  data. 
The  fact  that  the  experimental  data  does  not  appear  to  be 
exactly  approaching  the  origin  leads  one  to  believe  the 
deviation  may  be  due  to  experimental  error.  At  the  lower 
wicket  gate  angles,  the  probe  traverses  were  nearer  to  the 
trailing  edge  of  the  wicket  gates,  and  some  measurements  were, 
therefore,  made  with  the  probe  partially  in  the  wakes  of  the 
gates.  If  the  measurements  were  not  perfectly  centered  around 


the  wakes,  the  error  in  the  angle  measurements  would  not  be 
completely  canceled  out  by  the  averaging  process,  resulting  in 
one  possible  source  of  error.  An  additional  source  of  error 
could  be  the  angle  of  the  wicket  gates.  The  linkage,  which 
controls  the  angles  of  the  wicket  gates,  has  some  hysteresis, 
which  may  be  enough  to  allow  the  gates  to  be  at  a slightly 
different  angle  than  expected. 

Over  the  full  range  of  wicket  gate  angles  the  potential 
solution  appears  to  be  slightly  more  accurate  than  the  graphical 
method  for  the  symmetrical  wicket  gate  geometry.  However,  the  very 
small  difference  between  the  predictions  by  both  methods  and  the 


experimental  data  does  not  give  one  method  a significant  advantage 
over  the  other.  It  was  for  this  reason  that  an  analytical  study 
of  a system  with  cambered  wicket  gates  was  undertaken  to  investigate 
if  the  two  methods  differ  more  significantly  for  cambered  wicket 
gates.  The  geometry  of  the  wicket  gate  selected  for  the  study  is 
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pictured  In  Figure  18.  The  results  from  the  study  are  presented  In 
Figure  16,  where  the  fluid  angles  predicted  by  the  two  methods  are 
plotted  similarly  to  the  previous  data.  Figure  16  shows  a slightly 
larger  difference  between  the  two  methods  than  was  observed  for  the 
symmetrical  wicket  gate,  although  the  difference  is  still  not  extremely 
large.  For  the  cambered  wicket  gate  the  maximum  difference  is  approximately 
6 degrees  while  for  the  symmetrical  wicket  gate  the  difference  is 
approximatly  3.5  degrees.  For  the  cambered  wicket  gates  the  potential 
flow  solution  is  expected  to  yield  the  more  accurate  results,  since 
it  can  account  for  the  details  of  the  wicket  gate  shape.  The  significance 
of  this  conclusion  should  be  determined,  however,  and  cambered  wicket 
gates  should  be  studied  experimentally  in  the  future.  One  advantage 
of  the  potential  flow  solution  is  that  stay  vanes  can  be  added  to  the 
analysis  for  cases  where  the  gate  spacing  is  such  that  inlet  conditions 
to  the  gates  will  affect  the  exit  flow  angle.  The  pressure  distribution 
on  the  wicket  gates  can  also  be  obtained  from  the  potential  flow  solution, 
from  which  the  force  and  moment  coefficients  for  the  spindle  can  be 
calculated . 
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the  same  direction  as  the  mean  swirl  of  the  flow  and  opposite  to 
the  positive  6 direction  shown  in  Figure  12.  In  general.  Figure  17 
shows  that  the  predicted  and  measured  data  have  the  same  trends  and 
that  the  slope  of  data  is  almost  identical.  The  mean  value  of  the 
two  sets  of  data  is  the  main  difference,  which  tends  to  make  the 
prediction  look  worse  than  it  is.  The  mean  values  differ  only  by 
approximately  2 degrees.  It  should  be  pointed  out  that  the  large 
extreme  angles  in  the  experimental  data  are  erroneous  measurements 
made  with  the  probe  in  the  wake  of  wicket  gates.  Eliminating  these 
points  would  make  the  prediction  appear  more  accurate. 

3.8  Summary  and  Conclusions 

A two-dimensional  potential  flow  solution  has  been  presented 
for  the  flow  through  a radial  cascade,  such  as  the  wicket  gates  of 
a hydraulic  turbine.  Comparisons  of  the  predicted  flow  angles  from 
the  potential  solution  and  the  angles  measured  in  an  air  model  of 
a wicket  gate  system  indicate  that  the  potential  theory  very  adequately 
describes  the  real  fluid  characteristics.  For  the  wicket  gate  system 
investigated,  little  difference  was  found  between  the  potential  flow 
solution  and  a graphical  method  currently  used.  A purely  analytical 
study  of  cambered  wicket  gates  was  conducted  which  showed  a slightly 
larger  deviation  between  the  two  methods.  It  is  expected  that  the 
potential  solution  yields  the  more  accurate  results,  although  this 
conclusion  was  not  verified  experimentally.  For  the  present,  the  main 
advantages  of  the  potential  flow  solution  are  the  ability  to  obtain 
a more  detailed  solution  of  the  flow  and  the  ability  to  include 
upstream  effects  if  they  are  considered  to  be  influential.  The  pressure 
distribution  on  the  wicket  gates  is  also  obtained,  and  the  force  and 
moment  coefficients  for  the  wicket  gate  spindles  can  be  computed. 
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(Measured  from  the  radial  position) 

Figure  15  - Average  Flow  Angle  versus  Wicket  Gate  Angle  Obtained 
for  the  Symmetrical  Wicket  Gates  from  the  Potential 
Flow  Solution,  Graphical  Prediction  Method,  and 
Experimental  Data 
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Figure  17  - Local  Fluid  Angle  versus  Circumferential  Location 
Measured  Downstream  of  the  Wicket  Gates  and 
Predicted  by  the  Potential  Flow  Solution 
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Figure  18  - Schematic  of  a Wicket  Gate  and  the  Cascade  of  Wicket 
Gates  Describing  the  Program  Input  and  Output 
Parameters 
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Appendix  A 

Results  of  Preliminary  Tests  Indicating  Effect  of 
Test  Apparatus  Geometry  on  Surge  Pressure 

The  test  facility  for  the  draft  tube  surge  study  as  shown  in  Figure  3 
was  not  the  original  test  apparatus  considered  In  the  program.  The 
original  draft  tube  surge  test  facility  was  designed,  fabricated,  and 
assembled  as  shown  in  Figure  14.  The  main  difference  between  the  original 
test  facility  and  the  one  finally  used  for  testing  is  the  location  of  the 
air  supply  flow  fan.  In  the  original  test  facility,  (Figure  14),  air 
was  drawn  from  the  open  atmosphere  into  the  test  section  and  dumped  into 
a stilling  chamber  by  means  of  a flow  fan  downstream  of  the  draft  tube. 

In  the  final  test  facility,  (Figure  3),  air  was  blown  into  a stilling 
chamber  first,  then  passed  through  the  test  section  and  dumped  to 
atmosphere. 

Examination  of  the  sketches  of  the  two  facilities  would  seem  to 
indicate  that  the  important  flow  characteristics  (pressure  parameter  and 
frequency  parameter)  produced  in  either  test  facility  would  be  essentially 
the  same.  The  authors  felt  that  it  was  important  for  the  test  facility 
to  produce  relations  between  the  surge  parameters  and  the  momentum 
parameter  similar  to  those  obtained  by  Palde  [4].  Preliminary  tests 
in  the  original  test  apparatus  provided  a relation  between  frequency 
parameter  and  momentum  parameter  much  like  that  obtained  by  Palde. 

However,  pressure  parameter  as  a function  of  momentum  parameter  was 
essentially  constant,  which  radically  disagreed  with  Palde' s data. 

In  an  attempt  to  obtain  pressure  parameter  data  similar  to  Palde' s, 
various  pressure  pick-up  devices  were  tried  at  various  locations 
along  the  draft  tube.  After  repeated  failure  to  obtain  the  desired 
data  in  the  original  test  facility,  a test  facility  similar  to  that 
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used  by  Paldc  (i.e.,  pushing  air  through  the  draft  tube  and  dumping 
to  atmosphere)  was  constructed  as  shown  in  Figure  3.  Surging 
characteristics  obtained  in  the  modified  test  facility  agreed  very  well 
with  Palde's  data.  The  altered  test  facility  solved  the  problem  but 
raised  a question.  Why  was  the  pressure  parameter  much  lower  for  the 
original  test  facility  compared  to  that  obtained  in  Palde's  or  the 
final  test  facility?  Velocity  profiles  obtained  in  the  draft  tube  at 
L/D“4 .40  for  each  facility  at  various  momentum  parameters  were  compared. 
Figure  A.l  shows  that  the  velocity  profiles  are  quite  similar  for  both 
test  facilities. 

In  an  effort  to  answer  this  question,  a brief  series  of  tests  were 
performed  to  investigate  the  effects  on  pressure  parameter  of  dumping 
the  draft  tube  flow  into  a stilling  chamber  and  then  dumping  to 
atmosphere.  Tests  were  performed  on  both  the  cylindrical  and  elbow 
type  draft  tubes.  Tables  A1  and  A2  show  the  different  geometries  that 
were  investigated.  With  any  given  facility  geometry  tested,  a survey 
of  pressure  parameter  as  a function  of  momentum  parameter  was  obtained. 

The  flowrate  through  the  test  facility  was  measured  by  an  orifice 

meter  and  the  surge  pressure  was  measured  through  a dynamically  calibrated 

pressure  transducer  by  an  RMS  meter  as  earlier  described  in  this  report. 

Figure  A. 2 shows  the  results  of  this  investigation  for  the  cylindrical 
draft  tube  and  likewise.  Figure  A. 3 corresponds  to  the  Fontenelle  (elbow 
type)  draft  tube.  The  numbers  on  the  curves  in  Figures  A. 2 and  A. 3 
correspond  to  the  geometries  so  numbered  in  Tables  A1  and  A2.  Comparison 
of  curves  (1)  and  (2)  in  Figure  A. 2 shows  the  effects  of  a stilling 
chamber  below  the  draft  tube  on  pressure  parameter.  Looking  at  all 
the  curves  of  both  Figures  A. 2 and  A. 3,  it  is  obvious  that  the  use 
of  the  stilling  chamber  in  all  cases  reduced  the  pressure  parameter 
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Table  A.l 


1.  STILLING  CHAMBER  AND  AIR  SUPPLY  DOWNSTREAM  OF  DRAFT  TUBE 


4.  STILLING  CHAMBER  AND  AIR  SUPPLY  UPSTREAM  OF  DRAFT  TUBE 


DISCHARGE  TO 
ATMOSPHERE 


5. 


STILLING  CHAMBER  WITH  SCREEN  DOWNSTREAM  OF  DRAFT  TUBE.  UPSTREAM  AIR  SUPPLY 


DISCHARGE 
TO  ATMOSPHERE 
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: 


I.  NO  STILLING  CHAMBER.  UPSTREAM  AIR  SUPPLY 


2.  STILLING  CHAMBER  DOWNSTREAM  OF  DRAFT  TUBE.  UPSTREAM  AIR  SUPPLY 


ELBOW  DRAFT  TUBE 


3.  STILLING  CHAMBER  (SHORT  CONNECTION)  UPSTREAM  OF  DRAFT  TUBE.  UPSTREAM  AIR  SUPPLY 
„ SHORT  CONNECTION 


AIR 

SUPPLY 


DISCHARGE  TO  ATMOSPHERE 


4.  STILLING  CHAMBER  (LONG  CONNECTION)  UPSTREAM  OF  DRAFT  TUBE.  UPSTREAM  AIR  SUPPLY 
T ,LONG  CONNECTION 


AIR 

SUPPLY 


DISCHARGE  TO  ATMOSPHERE 
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magnitude.  Notice  that  screens  placed  across  the  stilling  chamber 
(curve  (5),  Figure  A. 2)  forced  the  pressure  parameter  relation  of 


the  original  test  facility  to  become  almost  identical  to  those  of 


Palde  and  curve  (2) 


Out  of  this  brief  study  comes  speculation  that  stilling  chambers 
placed  downstream  of  a turbine  runner  may  be  an  effective  means  of 
reducing  draft  tube  surge  in  hydroelectric  pump-turbines.  Further 
study  is  required  in  this  area. 


A second  complication  has  been  suggested  by  H.  Falvey  of  The 
Bureau  of  Reclamation.  This  involves  the  geometry  of  model  test 
loops  which  have  stilling  chambers  located  up  or  downstream  of  model 


pump-turbines.  The  water  level  and  geometry  of  these  chambers  could 


provide  surge  characteristics,  based  on  model  tests,  that  are 


inconsistent  with  those  of  the  prototype 
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^ORIGINAL  TEST 
FACILITY  (Fig.  14) 


MOMENTUM  PARAMETER, 


Figure  A. 2 


Pressure  versus  Momentum  Parameter  for  the  Cylindrical 
Draft  Tube  as  a Function  of  Test  Facility  Geometry 
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